How cold is cold dark matter? Small scales constraints from the flux power spectrum 

of the high-redshift Lyman-a forest 
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We present constraints on the mass of warm dark matter (WDM) particles derived from the 
Lyman-a flux power spectrum of 55 high- resolution HIRES spectra at 2.0 < z < 6.4. From the 
HIRES spectra, we obtain a lower limit of ttiwdm <^ 1-2 keV (2a) if the WDM consists of early 
decoupled thermal relics and ttiwdm ft 5.6 keV (2a) for sterile neutrinos. Adding the Sloan Digital 
Sky Survey Lyman-a flux power spectrum, we get mwDM ~ 4 keV and ttiwdm ~ 28 keV (2a) for 
thermal relics and sterile neutrinos. These results improve previous constraints by a factor two. 



Introduction. Warm dark matter (WDM) has been 
advocated in order to solve the apparent problems of 
standard cold dark matter (CDM) scenarios at small 
scales (e.g. most notably: the excess of the num- 

ber of galactic satellites, the cuspiness and high (phase 
space) density of galactic cores and the large number of 
galaxies filling voids. These and other problems could 
be alleviated if the dark matter (DM) is made of warm 
instead of cold particles. The main effect of the larger 
thermal velocities would be to suppress structures below 
the Mpc scale. However, poorly understood astrophys- 
ical processes governed by the baryonic component of 
galaxies, along with numerical, observational and theo- 
retical aspects [1, H, H[ > have also to be considered in order 
to reliably model the spatial distribution of dark matter 
at small scales. The Lyman-a absorption produced by 
the intervening neutral hydrogen in the spectra of dis- 
tant quasars (QSOs), the so called Lyman-a forest, is 
a powerful tool for constraining dark matter properties. 
It probes the matter power spectrum in the mildly non- 
linear regime over a large range of redshifts (z = 2 — 6) 
down to small scales (1 — 80 ft- -1 Mpc) 0. In previ- 
ous work, Ref. Q used two samples of high-resolution 
QSO Lyman-a forest spectra at z ~ 2.5 to set a lower 
limit of 550 eV for the mass of a thermal warm dark 
matter candidate (2 keV in case of a sterile neutrino). 
More recently, Ref. Q and Ref. [H, using the Sloan 
Digital Sky Survey (SDSS) QSO data set at higher red- 
shifts and different methods significantly improved this 
limit by a factor ~ 4. Among the possible WDM can- 
didates, the most promising appears to be a sterile neu- 
trino with a mass in the keV range, which could be part 
of many particle physics models with grand unification 
(e.g. [9(). Because of a non-zero mixing angle between 
active and sterile flavor states, X-ray flux observations 
can also constrain the abundance and decay rate of such 
DM particles (e.g. [T3] ) . The constraints from Lyman- 
a forest data and those from the X-ray fluxes of astro- 



physical objects together put considerable tension on the 
parameter space still allowed for a sterile neutrino parti- 
cle with the phase-space distribution proposed by Dodel- 
son & Widrow (DW) [IH, although other non-standard 
scenarios must be explored flil ]. Here, we will use a new 
large set of high-resolution Lyman-a forest spectra in or- 
der to improve limits on the mass of a WDM particle. 

Data sets. We use two different data sets: i) the 
high resolution HIRES data set presented in [3] which 
consists of 55 QSOs spanning the range 2.0 < z < 6.4; 
it) the SDSS Lyman-a forest data of McDonald et al. 
14j . which consists of 3035 quasar spectra at low res- 
olution (R ~ 2000) and low signal-to-noise spanning a 
wide range of redshifts (z = 2.2 — 4.2). We have calcu- 
lated the flux power spectrum from the HIRES data for 
4 redshift bins with median redshifts z = 2.5, 3.5, 4.5, 5.5 
and 20 logarithmically spaced bins in wavenumber span- 
ning 0.002 < k (s/km)< 0.287, Note, however, that we 
use only 12 bins in the range 0.003 < k (s/km)< 0.077 
for our present analysis (a total of 48 points). The to- 
tal redshift path of the HIRES sample is Az = 29.1. 
The redshift path of the lowest redshift bin is compa- 
rable to that of the LUQAS sample of high-resolution 
spectra (Az = 13) which we used previously [l5T |. while 
the other three redshift bins have smaller redshift paths 
of approximately equal length. Note that our previous 
work based on the LUQAS sample used only the range 
0.003 < k (s/km)< 0.03. For the SDSS data set we 
use 132 flux power spectrum measurements Pp(k, z) that 
span 11 redshift bins and 12 k— wavenumbers in the range 
0.00141 < k (s/km)< 0.01778 (roughly corresponding to 
scales of 5-50 comoving Mpc). Since the HIRES spec- 
tra have higher resolution than the SDSS spectra we can 
use the flux power spectrum obtained from the HIRES 
data to extend our analysis to smaller scales. We have 
removed the damped Lyman-a systems (DLAs) from the 
observed HIRES spectra but decided not to attempt to 
remove the metal lines. We rely on the measurement of 
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the contribution of metal absorption to the flux power 
spectrum at z = 2.13 by [IBj]. We use this estimate to 
correct our measurement of the flux power spectrum at 
z = 2.5, and apply a smaller correction for the measure- 
ment at z = 3.5. For the two highest redshift bins we do 
not apply any correction to our measurements of the flux 
power spectrum as the metal contribution should be very 
small at these redshifts. The contribution of the metal 
absorption to the flux power spectrum at the scales con- 
sidered here depends only weakly on wavenumber and is 
degenerate with the value of the assumed effective optical 
depth [l6[ . Marginalizing over the observed range of the 
effective optical depth in our parameter analysis will (im- 
plicitly) account for a possibly different contribution to 
the flux power spectrum than we have assumed. We have 
tested the effect of continuum fitting errors by calculating 
the flux power spectrum of each observed Lyman-a for- 
est spectrum at z > 4 for 100 continuum fit realizations 
where we have adjusted the continuum level of the ob- 
served spectra randomly by a factor 1+e with e in the 
range [—0.1,0.1]. This should be a reasonable estimate 
of the effect of continuum fitting errors at high redshift. 
At smaller redshifts the continuum fitting errors should 
be below 4%. The variations of the flux power for these 
100 realizations lies well within the statistical errors of 
the flux power spectrum of our full sample and estimates 
differ by less than 5% from the case with e = 0. We have 
decided not to try to explicitly account for continuum 
fitting errors. The covariance matrix of the flux power 
spectrum for the HIRES data set was calculated with a 
jack-knife estimator. 

Method. Modeling the flux power spectrum of the 
Lyman-a forest for given cosmological parameters down 
to the required small scales is not straightforward and 
accurate numerical supercomputer simulations are re- 
quired. Here, we model the flux power spectrum with full 
hydro-dynamical simulations using a second order Tay- 
lor expansion around a best fitting model. This allows 
us to obtain a reasonably accurate prediction of the flux 
power spectrum for a large range of parameters, based on 
a moderate number of simulations [17| . The method has 
been first introduced in Ref. [l6| and we refer to this work 
for further details. Note that while in Ref. [l6| the pre- 
diction for the flux power were made using a first order 
Taylor expansion, here the expansion is made to second 
order: i.e. the parameter dependence of the flux power 
spectrum Pp(k, z,p) is locally described by a 2nd order 
polynomial function for any redshift z, set of wavenum- 
bers k and cosmological or astrophysical parameters p. 
In this work, we use as input in the hydro-dynamical 
simulations a linear matter power spectrum as in [f| [25| 
and we assume that the sterile neutrino phase-space dis- 
tribution is equal to that of active neutrinos multiplied 
by a suppression factor. Deviations from this first-order 
approximation were computed in [l8| , but typically these 
corrections lower the bounds on the warm dark matter 



mass by only 10% 0]. 

For our best estimate of the flux power spectrum of 
our fiducial model we used a simulation of a box of length 
60 h^ 1 comoving Mpc with 2 x 400 3 gas and cold DM 
particles (gravitational softening 2.5 hr 1 kpc). The fidu- 
cial flux power spectrum has been corrected for box size 
and resolution effects. Note that resolution corrections 
are large at high redshift and redshift /scale dependent: 
they reach 50% (300%) at the smallest scales at z = 4.5 
(z = 5.5). We performed a number of additional hydro- 
dynamical simulations with a box size of 20 hr 1 comoving 
Mpc and 2 x 256 3 gas and DM particles (grav. soft. 1 
h^ 1 kpc) for WDM models with a (thermal) warm dark 
matter of mass wwdm = 1,4,8 keV, to calculate the 
derivatives of the flux power spectrum with respect to 
changes of the WDM particle mass and other astrophys- 
ical and cosmological parameters of interest in this anal- 
ysis. We have checked the convergence of the flux power 
spectrum on the relevant scales using several additional 
simulations (following the approach of [IH ) with 2 x 256 3 
gas and DM particles and box sizes of 10 h" 1 Mpc (grav. 
soft. 2 hr 1 kpc), 5 hr 1 Mpc (grav. soft. 1 h~ 1 kpc) and 
a 5 hr 1 Mpc simulation with 2 x 448 3 (grav. soft. 0.27 
h- 1 kpc) . 
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FIG. 1: Flux power spectrum of the HIRES data set at differ- 
ent redshifts and best fit models (solid curve) with biwdm = 8 
keV and a model with uiwdm = 2.5 keV (dashed curve). 

We then used a modified version of the code CosmoMC 
[20I ] to derive parameter likelihoods from the the HIRES 
and SDSS Lyman-a data. For the HIRES data, we had 
a set of 18 parameters: 7 cosmological parameters; 6 pa- 
rameters describing the thermal state of the Intergalac- 
tic Medium: parameterization of the gas temperature- 
gas density relation T = Tq(z)(1 + S) 1 ^" 1 as a broken 
power law at z — 3 with the two astrophysical param- 



3 



eters Tq (z) and j A {z) describing the amplitude; 4 pa- 
rameters describing the evolution of the effective optical 
depth with redshift (slope rjj and amplitude at z = 3 
and z = 5) since a single power-law has been shown to be 
a poor approximation over this wide redshift range (see 
[l3|); one parameter describing the spatial fluctuations 
of the Ultra- Violet (UV) background /uv- 

In estimating the effect of UV fluctuations on the 
flux power we adopt a conservative approach and con- 
sider the model of [2l|, where the UV fluctuation have 
a large impact on the flux. The model assumes that 
the UV background and its spatial fluctuations are pro- 
duced by Lyman-Brcak galaxies and QSOs and uses as 
input the QSOs and Lyman-Break luminosity functions 
at z = 3.5,4,5,6. At z = 5.5 the flux power in the 
model with UV fluctuations is larger by 4% at the largest 
scales increasing to 20% at k = 0.2 s/km (not included 
in the analysis), compared to the case without UV fluc- 
tuations. At z = 4 and z = 3.5 the only differences 
arise at scales k > 0.3 s/km, which are not considered in 
the present analysis. Further details on the UV model 
can be found in [2l| (but see other approaches [22I]). 
We parameterize the effect of UV fluctuations on the 
flux power with a multiplicative factor fuv constrained 
to be in the range [0,1]. For the SDSS data we have 
used a total of 28 parameters: 15 parameters used for 
the HIRES spectra (without /uv and the two param- 
eters describing the effective optical depth evolution at 
3 = 5) plus 13 noise-related parameters: 1 parameter 
which accounts for the contribution of DLAs and 12 pa- 
rameters modelling the resolution and the noise prop- 
erties of the SDSS data set (see [23|). We do not ad- 
dress the role of different reionization scenarios on the 
flux power. To do this self-consistently would require 
radiative transfer simulations beyond present numeri- 
cal capabilities and the effect of the reionization history 
should be subdominant and degenerate with the ther- 
mal state of the gas. In computing the likelihood a 
crucial input is the covariance matrix of the two data 
sets. The covariance matrix of the SDSS flux power is 
provided by the authors of 14j . We found the covari- 
ance matrix of our HIRES data set to be rather noisy 
(especially at high redshift), preventing a reliable inver- 
sion. To overcome this problem we use the suggestion of 
[24| . We regularize the observed covariance matrix using 
the correlation coefficients as estimated from the sim- 
ulated spectra, cov d (i,j) = r s (i,j)y/covd(i ,i)covd(J,j) 
withr s (i, j) = cov s cov s (i, i)cov s (j, j ) , where cov s 

and covd are the covariance matrices of the observed and 
simulated spectra, respectively. Note that this procedure 
implicitly assumes that observed and simulated data have 
similar covariance properties. We have applied moderate 
priors to the thermal history to mimic the observed ther- 
mal evolution as in [l|| and a prior on the Hubble con- 
stant (72 ±8 km/s/Mpc), but note that the final results 
for the mass constraint are not affected by these priors. 



Results. In Figurc[T]we show the best fit model for the 
HIRES data set (continuous curve, towdm = 8 keV) and 
a model with a smaller mass for the thermal WDM par- 
ticle (dashed line, towdm = 2.5 keV ). The constraining 
power of the small scales at high redshift is immediately 
evident. The % 2 value of the best fit model is ~ 40 for 36 
d.o.f. and with a probability of 16% this is a reasonable 
fit. As noted in Ref. 0] at high redshifts, the mean flux 
level is lower and the flux power spectrum is closer to 
the linear prediction making the flux power data points 
very sensitive to the free-streaming effect of WDM. We 
confirm that there are no strong degeneracies between 
mwDM and the other parameters, demonstrating that 
the effect of a WDM particle on the Lyman-a flux power 
is unique, and that the other cosmological and astrophys- 
ical parameters considered here cannot mimic its effect. 

The 2a lower limits for the mass of the warm dark 
matter particle are: 1.2 keV, 2.3 keV and 4 keV, for the 
HIRES, SDSS and SDSS+HIRES data sets, respectively. 
The corresponding limits for DW sterile neutrino are: 
5.6, 13, and 28 keV (see [f| for how the masses are related 
for the two cases). The x 2 °f the best fit model of the 
joint analysis ~ 198 for 170 d.o.f. which should occur in 
7% of the cases. The sample of HIRES spectra improves 
our previous constraint from high-resolution spectra ob- 
tained from the LUQAS sample by a factor two. Drop- 
ping the highest redshift bin [z = 5.5) weakens the limit 
to 0.8 keV (3.3 keV) for the mass of a thermal (sterile) 
neutrino. The SDSS data alone is still more constraining 
than the HIRES data alone, due to the smaller statisti- 
cal errors of the SDSS flux power spectrum and the finer 
coverage of a large redshift range which helps to break 
some of the degeneracies between astrophysical and cos- 
mological parameters. Combining the SDSS data and 
the HIRES results in an overall improvement of a factor 
~ 2 and gives the strongest limits on the mass of WDM 
particles from Lyman-a forest data to date. In Table 1 
we summarize the constraints obtained for the most rele- 
vant astrophysical and cosmological parameters (la) for 
our analysis of the HIRES only and HIRES+SDSS data 
sets. We note that, similarly to [f|, there is a preference 
for a non-zero 1 /?tiwdm value which is at present not sta- 
tistically significant (less than 2a). The data also prefers 
models with non-zero UV background fluctuations. 

Discussion. We have used a sample of high resolution 
Lyman-a forest spectra which is sensitive to the suppres- 
sion in the matter power spectrum at small scales caused 
by the free-streaming of WDM particles. We have mod- 
elled the observed flux power spectrum by using high 
resolution hydro-dynamical simulations that incorporate 
the relevant physical processes. We also improved previ- 
ous analyses by extending the parameter space, perform- 
ing a Taylor expansion order of the flux power spectrum 
in the cosmological and astrophysical parameters to sec- 
ond instead of first order, and including UV fluctuations 
that should be important at the high redshifts consid- 



4 



TABLE I: Mar$ 
parameter 



;inalized estimates 
HIRES+SDSS 



(la errors) 
HIRES 
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0.28 ± 0.09 
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3.17 ±0.07 


3.02 ±0.37 




= 3) 


1.44 ±0.12 


1.54 ±0.33 




= 5) 


1.53 ±0.09 


1.54 ±0.19 


-riff (z 


= 5) 


4.77 ±0.44 


4.92 ±0.5 


T (z-- 


= 3)(10 4 ) K 


2.23 ±0.30 


1.54 ±0.34 


/uv 




0.65 ±0.25 


0.58 ±0.28 


1/mwDM (keV 1 ) 


0.09 ±0.07 


0.44 ±0.22 



ered here. We confirm that the observed Lyman-a forest 
flux power spectrum at small scales and high rcdshifts 
requires significantly more power on small scales than 
provided by the models that try to reproduce the cores 
of dwarf galaxies with a warm dark matter particle by Q • 
We improve previous lower limits on the mass of warm 
dark particles of 0] by a factor two and those of Q by a 
factor three. This further decreases the rather small gap 
between the limits on the mass of DW sterile neutrinos 
from Lyman-a forest data and those on mass and mixing 
angle from the diffuse X-ray background (e.g. [III). We 
note that all the limits quoted here refers to the thermal 
production mechanism for sterile neutrinos and could be 
potentially 10-20% weaker in case of non-thermality (2(|. 
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